Draft version June 4, 2009 

Preprint typeset using I^T^X style cmulatcapj v. 10/09/06 



CONSTRAINING THE POROSITIES OF INTERSTELLAR DUST GRAINS 

Kevin Heng 1 & Bruce T. Draine 2 

Draft version June 4> 2009 

ABSTRACT 

We present theoretical calculations of the X-ray scattering properties of porous grain aggregates with 
olivine monomers. The small and large angle scattering properties of these aggregates are governed 
by the global structure and substructure of the grain, respectively. We construct two diagnostics, IZx 
and 7x, based on the optical and X-ray properties of the aggregates, and apply them to a Chandra 
measurement of the dust halo around the Galactic binary GX13+1. Grain aggregates with porosities 
V > 0.55 are ruled out. Future high-precision observations of X-ray dust haloes together with detailed 
modeling of the X-ray scattering properties of porous grain mixtures will further constrain the presence 
of porous grain aggregates in a given dust population. 

Subject headings: ISM: dust, extinction — scattering — X-rays: diffuse background 



1. INTRODUCTION 

After many decades of trying to deduce the properties of interstellar dust from observations, the basic geometric 
structure of the grains remains controversial. The observed polarization of starlight establishes that some or all grains 
are non-spherical, and there is agreement that interstellar grains must include both silicate material and carbonaceous 
material, but there are divergent views concerning the internal structure of grains. Some authors (e.g., Mathis et al. 
1977; Draine & Lee 1984; Kim, Martin & Hendry 1994; Kim & Martin 1995; Weingartner & Draine 2001, hereafter 
WD01; Draine & Li 2007; Draine & Fraisse 2009) have modeled interstellar grains as compact, zero-porosity spheres 
or spheroids, with some consisting of silicate and others consisting of carbonaceous material. Others have argued that 
the evolution of interstellar grains involves coagulation in clouds, producing grains with a mixed composition and a 
"fluffy" structure, with porosities taken to be ~ 0.8 (Mathis & Whiffen 1989) or > 0.9 (Voshchinnikov et al. 2006). 

Because polarized light-scattering is sensitive to the geometry of the scatterers, it has been possible to use the 
observed color and polarization of scattered light, as a function of scattering angle, to estimate the size and porosity 
of cometary dust particles and the dust in debris disks. Shen et al. (2008, hereafter SDJ08) examined the optical 
properties of porous grain aggregates using the discrete dipole approximation (Draine 1988; Draine & Flatau 1994). 
Shen et al. (2009) found that moderate-porosity aggregates can reproduce the observed scattering properties of dust in 
the AU Mic debris disk and in comets. Unfortunately, the scattering properties of interstellar grains remain uncertain, 
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Fig. 1. — Random aggregates produced by the BAM2, BAM1, and BA agglomeration schemes (Shen, Draine & Johnson 2008). The unit 
vector Si is along the principal axis of the largest moment of inertia. With ag = 0.04/im silicate spherules, each of these aggregates (with 
a cS = 0.254, 0.320, and 0.508/mi) has Q e xt (-B)/Q cx t (R) ~ 1.70, the observed value of A B /A R (see text). 
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Fig. 2. — Plots of dQ BC!i t /dSl (E = 2 keV) for (left) a single spherical grain with radius a = 0.25 fim; and (right) the BAM2 aggregate of 
Fig.[T] with 256 ao = 0.04 fim silicate spherules (P = 0.563, a e g 0.254 (im), oriented with 9 = 0°, after "/3-averaging" (see text). 

and to date there has been no observational diagnostic capable of discriminating between these different models for 
interstellar grains. 

The present study shows that X-ray scattering can discriminate between compact and fluffy models for interstellar 
dust. Compact grain models have been successful in reproducing observations of X-ray scattering halos (e.g., Draine & 
Tan 2003, Smith 2008). In the present paper, we show how the X-ray scattering properties of irregular porous grains 
can be accurately calculated. We find that compact and fluffy grains differ substantially in their X-ray scattering 
properties, because the X-rays can "see" the small-scale structure within an aggregate grain, and because fluffy grains 
must be larger than compact grains when required to reproduce observed interstellar reddening constraints. 

In fJ2J we describe the X-ray optics of grain aggregates, including approximate scaling laws and our application 
of anomalous diffraction theory (ADT) to calculate their X-ray scattering properties. We present the results of our 
calculations in fJU In fj4j we construct two diagnostics based on the optical and X-ray properties of the dust grains and 
apply them to the Chandra halo measurement around the Galactic binary GX13+1 by Smith (2008, hereafter S08). 
We are able to rule out grain aggregate models with porosities "P>0.55. We summarize and discuss our results in Sj5] 

2. X-RAY OPTICS OF GRAIN AGGREGATES 

2.1. Anomalous Diffraction Theory (ADT) 

ADT, a combination of ray-tracing optics and Huygens' principle of propagation (van de Hulst 1957), is applicable 
to the study of dust grains at X-ray energies, when the refractive index is nearly unity and the grain size is relatively 
large (van de Hulst 1957): 

|to-1|«1, 
k a > 1, 

where m is the complex refractive index, Ao = 27r/fco is the wavelength of the incident photon and a is the size of 
the dust grain. Draine & Allaf-Akbari (2006, hereafter DA06) applied the ADT to calculate X-ray scattering from 
non-spherical grains. The interested reader is referred to DA06 or chapter 10 of van de Hulst (1957) for more details. 

Let ki = (0,0, fco) and k s = (k x , k y , ^Jk^ — k% — k%) be the incident and scattered propagation vectors. The differ- 
ential scattering cross section is given by 

where the scattering function, 

S (k x ,k y ;k ) = k„ T{f (x,y;k Q )} , (3) 
is related to the Fourier transform of the shadow function, / = f(x,y; k ): 

F {/} = J exp[i(k x x + k y y)} f dxdy. (4) 

The shadow function quantifies the fractional change in the propagating electric field, approximated as a plane wave, 
at a plane located just beyond the grain: 



f(x,y,k ) = 1-exp 



ifc / (to — 1) dz 



(5) 



where m = m(x, y, z) in general. 
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2.2. Scaling Laws 

We define the "effective radius" (or "volume equivalent radius"), 

3^ 1/3 



a cff =^j , (6) 

where V is the volume of solid material in the grain. If V is the grain porosity as defined by SDJ08, then the 
characteristic size is 

R = a cB (l - V)- 1 / 3 . (7) 

1 /3 

The grain is assumed to be made up of N s spherical monomers ("spherules") with radii do such that a e s = N s ao- 
The porosity V = V(N S ), depends on the agglomeration scheme used (see §3}- In Fig. [TJ we show visualizations of 
three random aggregates. 

Consider an X-ray photon with wave number ko encountering a dust grain and deflected by an angle Q and a 
transverse wave number, k± = ko sinO. The characteristic scattering angle is 

6char = JL_ _ 

k R irR 

(°) 



271" {1-V) 1/3 



1 keV\ f 0.3pm 



E J \ a c fr 



The total scattering cross section is 

O-scat « y—^f-J _ ^char- (9) 

We can derive approximate scaling laws for da sca t / dn at small and large angles. At sufficiently high energies such 
that |/| -C 1, we have 

•H/}«E ex p(^) ^ ( 10 ) 

3 

where To is the Fourier transform of the shadow function for one spherule and <pj is the phase contribution by the j-th 
spherule. 

At small scattering angles (k±R < 1; "core"), we have To « —ikQ(m — 1)Vq where Vq — Aira^/i is the volume of one 
spherule. The phase shifts are small {{4> 2 ) < 1) and coherent scattering occurs: 

T « N S T . (11) 

It follows that 

^«(AWo 2 )V-l| 2 . (12) 
At large scattering angles (k±R ^> 1; "wing"), the phase shifts are large and pseudo-random so that 

_ \T\*n¥ 2 \T \. (13) 
To evaluate To, we examine the scattering function for a single spherule (van de Hulst 1957): 

2 r' 2 

So (6; ko) = (fcoao) / [1 — CX P ( — ipsinu)] Jo (x cos it) sinw cos u du, (14) 



(15) 



where p = 2koao(m — 1), x = fco a o© an d Jo is the zeroth order Bessel function. When x ^> 1 and p ~ 1, 

So (6; ko) w * (|) V2 (fcoao) 2 PX" 3/2 ^3/2 (x) 



since (Arfken & Weber 1995) 
for x 3> 1. We then have 



2? (m — 1) fcoflo© 2 

x- ^ ( - 



«7i (x) ~ \ — cos 

V 7TX 



(16) 



da* 



/n «4iV s |m-l| 2 a2e- 4 = 4|m-l| 2 p£ 6- 4 . (17) 
ail \ ao J 

The scattering at large angles is proportional to the "power" at large k± in \T(f)\ 2 . 
Consequently, the normalized differential scattering cross section follows the scaling laws: 

rfQscat _ 1 da scat \Ns /3 ,0«Ochar, /.gs 

dn ~ ira 2 cS dn K j^e- 4 ,e>e char . 

We expect these scaling laws to be obeyed when absorption effects are weak: 

r = koa cSmi (1 - P) 2/3 < 1, (19) 
where is the imaginary part of the index of refraction. For E > 1 keV, we have m, < 1.9 x 10~ 4 for olivine 
(MgFeSi0 4 ) and thus r < 0.3(1 - V) 2/3 for a cS = 0.3 pm. 
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Fig. 3.— (eiQsart/dfi) at E = 2keV as a function of the scattering angle, 6, for BA.256.1 (V = 0.860), BAM1. 256.1 (V = 0.706) and 
BAM2. 256.1 (V = 0.563) aggregates, each consisting of 256 ag = 0.04 /an silicate spherules. 
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Fig. 4.— Comparing ADT and EMT calculations for (dQ scat /dU) at E = 2keV for BAM2.256.1 with a = 0.04 /im. 

2.3. Implementation 

To calculate the differential scattering cross section in equation (j2j, we adapt the Fortran 90 routine four2.f90 
from Press et al. (1996) for computing two-dimensional fast Fourier transforms (FFTs) with complex input functions. 
We are able to reproduce all of the results in Fig. 1 of DA06, thus verifying the accuracy of our code. As in the case 
of DA06, we find that a square, cartesian grid (x, y) with (2 11 ) 2 to (2 12 ) 2 points suffices for our purposes and avoids 
spurious contributions to the FFT resulting from aliasing. Bilinear interpolation is used to map the calculations into 
polar coordinates. 

3. RESULTS 

SDJ08 studied three types of grain aggregate models with different ranges of porosities. The most porous grains are 
the BA ("ballistic agglomeration") aggregates, constructed by requiring the arriving monomers to adhere to the point 
where they first made contact with the pre-existing aggregate. Less porous are the BAM1 ("BA with one migration") 
aggregates, where an arriving monomer is required to roll or slide along the contacted monomer, via the shortest 
possible trajectory, until a second contact point is established with the aggregate. The most compact aggregates 
are the BAM2 ("BA with two migrations"). For 256 < N s < 4096, the BA, BAM1 and BAM2 aggregates have 
(V) ~ 0.85-0.86, 0.74-0.78 and 0.58-0.66, respectively (see Table 2 of SDJ08). The current study employs random 
realizations of BA, BAM1, and BAM2 aggregates. 3 The monomers are assumed to have composition MgFeSi04 and 
density p — 3.8 gem - of olivine, with refractive index from Draine (2003). 

The orientation of each aggregate is defined by two angles: 9 is the angle between the line of sight and the principal 
axis of the largest moment of inertia, with unit vector <X\) and (3 is the rotation angle about a\. Rotation of a\ about 

3 Available in electronic form at http://www.astro.princeton.edu/~draine/agglom.html. 
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Fig. 5. — 7^x vs. the ratio of B- to ij-band extinction cross sections. Higher 7?.x and lower <jb/ijr values correspond to larger grain 
sizes. We have considered grain aggregates consisting of N s = 256, 512, 1024, 2048 and 4096 silicate spherules with radii ao = 0.04 (im. 
For comparison, calculations for single silicate spheres are shown. Also shown is the X-ray halo calculated for the WD01 mixture of silicate 
and carbonaceous grains. The dotted vertical line is crg/crR = 1.689 observed for typical diffuse cloud sightlines with Ay/E(B — V) fa 3.1 
(Fitzpatrick 1999). 

the line of sight can be suppressed because the scattered intensity will later be azimuthally averaged. Because of the 
reflection symmetry of the shadow function, we need to calculate S only for < 9 < ir/2 and < (3 < n. For 11 values 
of 9, we average dQ sca t/d£l over 11 values of (3 ("/?- averaging"). 

We compute dQ sca t/dtt for the BAM2. 256.1 aggregate with N s = 256 in the right panel of Fig. [2] The effective 
radius for this grain aggregate is a g = 0.254 fim (i? w 0.335 fim). For comparison, we calculate dQ scgut /dD, for a solid 
sphere with a — 0.25 fim. The first feature to notice is the series of distinct rings in the case of the solid sphere. These 
are the two-dimensional analogue of the Fourier transform of a Heaviside function in one dimension, and are artifacts 
of the sharply-defined edge of the spherical grain. For the aggregate, traces of these rings persist, but they are smeared 
out because there is no sharp spherical outer edge. 

Secondly, the scattering has similar intensities for both the spherical grain and the grain aggregate, but the central 
peak is concentrated within a smaller angular area for the latter. This is a simple consequence of the fact that porous 
grains generally have larger characteristic sizes than solid spheres. Larger grain sizes result in more concentrated 
forward scattering. 

For the rest of the paper, we average dQ sca t/dQ, over 11 values of the orientation angle 9 ("^-averaging"). We also 
define (dQ sca t/ d£l) to be the azimuthal average of the two-dimensional dQ scai / dVL function. We show examples of 
(dQscnt/dfl) for N s = 256 aggregates with BA, BAM1 and BAM2 porosities in Fig. [3j As expected from equation (fl"2|) . 
the BA, BAM1, and BAM2 aggregates have nearly identical forward scattering, but the BA clusters, being "larger", 
have a narrower forward scattering lobe, and therefore weaker scattering at intermediate angles (80-200" in Fig. [3]). At 
large angles 9 > 500" the BA, BAM1, and BAM2 clusters in Fig. [3] (all composed of the same number N s of identical 
spherules) have similar da sca ±/dQ. For grains with different values of N s , we checked the scaling laws described in 
equation (fT8|) and verified that they are obeyed to within a few percent. 

In Fig. 21 we include as an example the effective medium theory (EMT) calculation for a grain with the same mass 
and porosity as a N s = 256, BAM2 aggregate. The EMT calculation approximates the cluster by a uniform density 
sphere with refractive index tuemt = 1 + (1 — V){m — 1). For 9 < O c har, the EMT and ADT calculations are in 
agreement. However, EMT does not take into account the substructure of the grain aggregate, which is important for 

e > e char . 

We have also used the discrete dipole approximation (Draine & Flatau 1994) to calculate the extinction cross sections 
for the clusters at the effective wavelengths of the R, V, and B bands (A = 0.6492, 0.5470, and 0.4405 /im), using the 
public-domain code DDSCAT 7.0 (Draine & Flatau 2008). Results for selected clusters are shown in Table [TJ together 
with the differential cross section for forward-scattering of X-rays. 

4. X-RAY SCATTERING DIAGNOSTICS 

The scattering angle 8 is related to the observed halo angle ^ via the relation: 

d tan * = (D - d) tan (9 - (20) 

where d and D are the distances to the dust population and X-ray source, respectively. For example, we have 'J = 9/2 
when d — D/2. Generally, we have \l/ when d/D <C 1 and $ < 8 when d ~ D. 

For a e fj>0.1 /xm grains at optical wavelengths, the ratio obI&r of extinction in the B- to the i?-band is a generally 
decreasing function of a c g and serves as a diagnostic of the grain size. The observed value of Ab/Ar w 1.689 for 
the average interstellar reddening curve with Ay/E(B — V) ks 3.1 (Fitzpatrick 1999) requires viable grain models to 
have the extinction at optical wavelengths dominated by grains with (Jb/^b. ~ 1.689. Single-sized spheres with the 
dielectric function of "astronomical silicate" (Draine 2003) have <jbI<Jr = 1.689 for radius a = 0.177^m. 
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Fig. 6.— [e 2 dQ sca (2keV)/dQ]/Q ext (y band) for BA. 2048. 1-3, BAM1.512.1-3 and BAM2.256.1-3 aggregates (with a = 0.04 /im) that 
approximately satisfy crg/an = 1.689 (see text). Forward scattering per unit Ay increases with increasing porosity. The result for the 
WD01 dust model is also shown. 

A grain model can be tested by comparison with observed X-ray scattering halos over a range of halo angles. Because 
calculations of optical-UV extinction and X-ray scattering by complex grain geometries are very time-consuming, it is 
helpful to have simple diagnostics that can be applied specifically to the grains that dominate the extinction at optical 
wavelengths, because these same grains are expected to dominate the total X-ray scattering at all except the largest 
angles. Accordingly, we seek to characterize the X-ray scattering properties of the grain model for scattering angles 
e<e char (a = 0.2/rai) « 407" (keV/E). 

4.1. Core Slope Diagnostic: 72x 

For interstellar dust models based on compact grains, the optical extinction is dominated by grains with radii 
R « 0.2/im, which will have characteristic scattering angles c har ~ 200"(2 keV/E). For dust at d/D « 0.5, the 
characteristic halo angle would be \& c har ~ O.50 c har- Here we devise a diagnostic 72.x that is based on the slope of the 
scattering halo for \I> < Cellar: 

72 x (*!, * a , * 3 , * 4 ; E U E 2 ) = ^^l^l 1 ] , (21) 

Act (W3, W 4 , £j2) 



where 



ACT(* a ,vI/ b ,i?) = 27r F" *?£*(Q,E) sin*rf*. (22) 

dil 



For fixed halo angles \&i < ^2 < ^3 < ^4, we have higher values of TZx for larger grains because of the increased 
forward scattering. An integral property like 72x provides a useful comparison between theory and observations because 
it averages over oscillatory behavior in e?CT sca t / dCl(Q) when considering grains of a single size. 
For illustration, we set E\ = E<z = 2 keV, and adopt the following halo angles: 

# 1 = 5" (< e char ) 

* 2 = * 3 = 50" (» 0.259 char ) , (23) 
* 4 = 100" (« 0.58 char ) 

In Fig. El we construct curves of IZx vs. cfb/pr, f° r different values of d/D. The grain aggregates considered have 
0.25 < a s < 0.64 /im and 0.33 < R < 1.25 /im. For the calculations in Fig. we have averaged the differential cross 
sections over three random realizations of the grain aggregates, a process we will incorporate for the rest of the paper. 

The vertical dotted line in Fig.[5]is the "observed" value of <jbI&r = 1-689 for dust in the diffuse interstellar medium 
(ISM) with Ry = Av/E(B — V) = 3.1 (Fitzpatrick 1999). It represents an average over the mixture of grain sizes 
and composition present in the ISM. For single-size grains with a given porosity, ctb/ctjj decreases monotonically with 
increasing size for ctb/ctr > 1.2. The interstellar grain size distribution extends over a wide range, but the grain size 
for which (Jb/ctr = 1.689 will be characteristic of the grains that dominate the extinction in the B-, V- and i?-bands. 

The computations reported for the aggregates were CPU-intensive, therefore only a limited number of cases could be 
calculated. In Fig.[6l we show the differential scattering cross section for the BA (N s — 2048), BAM1 (N s = 512) and 
BAM2 (N s = 256) grains that approximately satisfy the ctb/ctr = 1.689 constraint. We compare these calculations to 
the scattering properties of the WD01 dust model. 

72.x can be calculated for any grain model. However, grain models that are incompatible with the observed interstellar 
reddening law are of no interest. Restricting ourselves to aggregates that satisfy ctb/cta = 1.689 determines N s for our 
chosen value of ao = 0.04 /im. 
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Fig. 7. — Observed and theoretical 7?-x vs. the ratio of the distances to the dust population and GX13+1. The width of the "observed" 
band reflects uncertainties in using both the HRCT and ACIST data to evaluate Atr (50", 100", E). 

We compare our results in Fig. [5] to calculations for single astronomical silicate spheres with radii of a = 0.1-0.3 
fim. The characteristic grain radius that satisfies (tb/cr, = 1.689 is a = 0.177 /xm. We also show IZx calculated for the 
WD01 size distributions for carbonaceous (including polycyclic aromatic hydrocarbons) and silicate grain populations 
for sightlines with Ry = 3.1. For d/D < 0.9, our single sphere calculations underestimate the IZx values calculated 
for the WD01 size distributions, though by only ~33% on average. 

Another issue to address is the grain composition adopted. We have neglected carbonaceous material, because the 
volume of silicate material is ~1.5 times the volume of carbonaceous solids. Moreover, the silicates are more effective 
for X-ray scattering at E > 1 keV because \m — 1| is larger than for graphite (see Figs. 3 and 5 of Draine [2003]). 
At E = 2 keV, \m — 1| = 1.94 x 10~ 4 for MgFeSiC>4 versus 1.13 x 10 -4 for graphite; the scattering power scales as 
\m — 1| 2 . Therefore, the graphite grains make only a secondary contribution to the scattering at E > 1 keV. 



We construct a second diagnostic based on the strength of the X-ray scattering in the "core" relative to the optical 
extinction. Define 



where Ar sca t = (nd/nH_)NnAa, {rid/nu) is the number of dust grains per H nucleon, Ah is the column density, and 
Act is defined in equation (|22[) . Tx can be calculated for any mixture of grains, but we are again only interested in 
grain models that are compatible with the observed interstellar reddening law. The forward-scattering cross section is 
strongly dependent on the grain size, and therefore this diagnostic is sensitive to the presence of larger grains in the 
size distribution. 



The IZx and Tx diagnostics constructed in § £14.114.21 allow us to compare our theory with observations of X-ray 
dust haloes. Recently, high-precision Chandra observations of scattered haloes around X-ray binaries have allowed for 
detailed comparisons to grain models (S08; Thompson & Rothschild 2009). S08 gives a flavor of the technical difficulty 
of the X-ray background subtraction involved in measuring the dust halo around the Galactic binary GX13+1 using 
HRC-I and ACIS-I. The HRC-I data is measured for 2" < * < 786", while the ACIS-I data is for 50" < * < 989". 
S08 noted the good level of agreement between the HRC-I and ACIS-I data between 50" and 100". Uncertainties 
with background subtraction make it difficult to study X-ray halos at large halo angles (\I/ > 500"). At small angles 
(\l/ < 5"), it becomes difficult to untangle the halo signal from the instrumental point spread function (resulting from, 
e.g., micro- roughness of the telescope mirror). 

The range of halo angles available in the GX13+1 measurement makes it an ideal case study for us. However, the 
HRC-I count rate is dominated by the point-spread function (PSF) for > 200" (see §3 of S08) and the measured PSF 
provided to us by R.K. Smith (2008, private communication) in electronic form only extends up to about 180". Since 
the HRC-I has no energy resolution, S08 used the best-fit model for the measured RXTE PCA spectrum of GX13+1 
and folded it through the HRC-I response function (see Fig. 3 of S08). We compute the effective Ac X I , 2, E) by 



4.2. Core X-Ray Scattering per Ay Diagnostic: Tx 




(24) 



4.3. Comparison to Measured Dust Haloes 
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Fig. 8. — Observed and theoretical T>c vs. the ratio of the distances to the dust population and GX13+1. The width of the "observed" 
band reflects uncertainties in Ay. 

taking an average over the energy spectrum of detected counts: 

4 

A^{* 1 ,* 2 ,E) = J2WE j Av(*i,*2,E j ), (25) 
4=1 

where the index j runs over Ej — 1.75 — 3.25 keV in intervals of AE = 0.5 keV. The characteristic energy is then 
E w 2.5 keV; we set *i = 4" and # 2 = 40". The weights are 

fEj + AE/2 
w JEj-AE/2 N E dK 

W * = J E N E dE > (26) 

where N E is the energy spectrum of detected counts. Similarly, we compute Act (^3, ^4, E) by taking ^3 = 50" and 
* 4 = 100". Using Act * 2 , E) and Act (* 3 , * 4 , E), we construct theoretical values of Tlx in Fig. [7] 

Using the HRCT and ACIS-I data and the stated halo angles yields the measured value of 7?x,obs = 1-26 ± 0.25. 
The theoretical value depends on the assumed value of d/D. From Fig. [3 we see that the BA model appears to be 
ruled out by the Tlx test. The BAM1 and BAM2 models are compatible with the Tlx test for d/D s=s 0.25 and 0.5, 
respectively. 

We show theoretical values of Tx in Fig. [5J To obtain Tx.obs, we need Ay. Ueda et al. (2004) measured the column 
density along the line of sight to GX13+1 to be Ah = (3.2 ± 0.2) x 10 22 cm" 2 . Since GX13+1 is located towards 
the inner galaxy, we assume Ay /Ah to be 1.25 ± 0.25 times the local value (4.65 x 10 _22 magcm 2 ; Rachford et al. 
2009) and obtain Ay = 18.6 ± 3.9, somewhat larger than earlier estimates Ay = 10.2 ± 0.3 (Garcia et al. 1992) and 
Ay = 15.4 ± 2.2 (Charles & Naylor 1992). The observed value Tx.obs = 0.0046 ± 0.0013 is inconsistent with the BA 
and BAM1 models for all values of d/D, but may be marginally compatible with the BAM2 model provided d/D<0.1. 
However, it is clear that the large extinction toward GX13+1 is not located at d/D 0, but instead almost certainly 
arises mainly in the 3.5 kpc ring at d ss 5 kpc where the Norma and Crux-Scutum arms cross the sightline to GX13+1. 
Therefore we conclude that all 3 fluffy-grain models - BA (V « 0.87), BAM1 (V w 0.73), and BAM2 (V w 0.57) - are 
ruled out by 7x,obs toward GX13+1. Based on both diagnostics, the porosity of interstellar dust is V < 0.55. 

What do the Tlx and Tx diagnostics say about compact grain models? The WDOl model, based on a mixture 
of compact silicate and carbonaceous grains, reproduces 7?x,obs if d/D « 0.60 — 0.65, but reproduces 7x,obs only if 
d/D sa 0.2-0.4. In reality, the dust is presumably distributed along the sightline to GX13+1 rather than residing in a 
single cloud. From the variation of Tx with d/D in Fig. [8jit is clear that WDOl dust uniformly-distributed along the 
sightline would result in a Tx value considerably larger than 7x,obs, consistent with the conclusion of Smith (2008) 
who found that smoothly-distributed WDOl model dust would approximately reproduce the observed X-ray scattering 
only if A H ~ (1.5 - 2) x 10 22 cm~ 3 , well below A H ~ 3.2 x 10 22 cm" 2 estimated by Ueda et al. (2004). This may 
indicate that the WDOl dust model produces too much X-ray scattering. However, (1) there may be systematic errors 
in the modeling introduced by the broad spectrum of X-ray energies contributing to the HRC-I imaging, and (2) the 
effects of multiple scattering are not negligible at the lowest energies. In addition, the estimate of Ay is based on X-ray 
absorption along the line of sight to GX13+1, whereas the X-ray scattering at 300" is produced by dust that may be 
displaced from the line of sight by up to ~10 pc. Perhaps Ay on the direct sightline to GX13+1 is larger than the 
average dust column contributing to the scattered halo. For example, if the dust producing the scattering halo toward 
GX13+1 had Ay ss 12 ±3 rather than 18.6 ±3.9, the "observed" band in Fig. [5] would shift upward by about a factor 
of 1.55, and the WDOl model would be more-or-less consistent with both 7?-x,obs and 7x,obs for d/D w 0.55. Compact 
grain models therefore appear to be viable, whereas all three of the fluffy grain models considered here appear to be 
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TABLE 1 

Optical Extinction and X-Ray Forward Scattering for Selected Targets 



target 


T> 

r 


<*eff 


D 


O * 

Vcxt 


Vcxt 


Vcxt 


(rid /rIO\n 






( n,xx\\ 


{ IIXU \ 


R 


v 


B 


2 keV 


BA. 2048.1 


0.8644 


0.508 


0.989 


4.986 


6.484 


8.352 


3.353 x 10 6 


BA.2048.2 


0.8657 


0.508 


0.992 


4.958 


6.460 


8.370 


3.384 x 10 6 


BA.2048.3 


0.8651 


0.508 


0.990 


5.004 


6.502 


8.392 


3.364 x 10 6 


BAM1. 512.1 


0.7239 


0.320 


0.491 


3.358 


4.425 


5.764 


5.319 x 10 s 


BAM1.512.2 


0.7530 


0.320 


0.510 


3.211 


4.308 


5.749 


5.326 x 10 s 


BAM1.512.3 


0.7286 


0.320 


0.494 


3.331 


4.398 


5.789 


5.279 x 10 5 


BAM2. 256.1 


0.5632 


0.254 


0.335 


2.877 


3.771 


4.877 


2.116 x 10 6 


BAM2.256.2 


0.5781 


0.254 


0.339 


2.841 


3.736 


4.822 


2.103 x 10 s 


BAM2.256.3 


0.5818 


0.254 


0.340 


2.807 


3.732 


4.857 


2.085 x 10 s 


sphere 





0.177 


0.177 


2.499 


2.940 


4.227 


5.012 x 10 4 



All targets consist of astronomical silicate (Draine 2003). 



firmly ruled out by the Tx diagnostic. 4 

5. DISCUSSION 

We have shown how the global structure of porous dust grain agglomerates leads to increased X-ray scattering at 
small angles, when compared with the scattering properties of spherical grains (both single size and the WD01 size 
distribution) that are constrained to have the observed ratio of extinction in the B and R bands. The enhanced small- 
angle scattering implies that observations of scattered X-ray halos can be used to constrain the internal geometry of 
interstellar grains. Applying the TZx and Tx diagnostics described above, we find that the observed X-ray halo around 
GX13+1 is inconsistent with grains consisting of random aggregates with porosities V > 0.55, at least for the monomer 
sizes considered here (ao = 0.04 ^m). 

Smith et al. (2002) argued that the dust toward GX13+1 could not consist of porous grains, on the grounds that 
porous grains could not produce sufficient overall scattering. Our Tx diagnostic uses the strength of the "core" of the 
scattered halo to constrain the dust porosity, but now we use the fact that porous grains would produce too much 
small-angle scattering per unit Ay. 

This exploratory study has focused on the X-ray scattering properties of only the grains that dominate the optical 
extinction, with the expectation that these are representative. Future studies should: 

1. Use size distributions of grain aggregates (adjusted to reproduce the full extinction curve and not just Ab/Ar) 
rather than just averaging over random realizations of a single aggregate size as was done here. We have shown 
that the TZx and Tx diagnostics for single-size silicate spheres with <tb/&r = 1-69 provide a good approximation 
to TZx and Tx calculated for the full WD01 size distribution, so we are confident that using single-size fluffy grains 
with ctb/(t_r ~ 1.69 provides a reasonable approximation to what TZx and Tx would be for size distributions 
constrained to reproduce the observed extinction curve. Nevertheless, TZx, m particular, will change when the 
single "representative" grain size is replaced by a size distribution that reproduces the IR-UV extinction curve. 
The Tx diagnostic is expected to be more robust, because the grains that dominate the visual extinction arc also 
expected to dominate the scattering at 6 < O c h ar (-R = 0.2/um). 

2. Include carbonaceous material, even though that material contributes only ^40% of the total grain volume. 
Again, we have already tested for sensitivity to the carbonaceous component by showing that TZx and Tx for 
single-size silicate spheres approximates TZx and Tx for the full WD01 size distribution, including carbonaceous 
grains, but it would be best to explicitly consider fluffy grain models that include a carbonaceous component. 

3. Examine the effect of varying the monomer size. The monomer size used in the present work (0.040 /im) was 
chosen arbitrarily. While we expect that similar results will be obtained for monomer sizes ao > 0.010 /um, this 
remains to be verified by future (numerically challenging) calculations. 

Although confirmation by more extensive future studies (as discussed above) will be invaluable, the diagnostics 
developed here, in particular the Tx diagnostic, already demonstrate that the dust toward GX13+1 has porosity 
"P<0.55. This finding appears to rule out models (e.g., Mathis & Whiffen 1989; Voshchinnikov et al. 2006) in which 
the interstellar grain mixture is dominated by highly-porous aggregates with V > 0.8. 

The porosity of interstellar dust grains therefore appears to be lower than the porosity of the submicron grains in 
the debris disk around AU Mic, which (based on interpretation of observed polarized scattering) have V rts 0.6, or 
micron-sized cometary dust grains, which also appear to have V w 0.6 (Shcn et al. 2009, and references therein). 

Why interstellar grains are not highly-porous is an open question. Gyroresonant acceleration by MHD turbulence in 
diffuse clouds appears able to drive a > 0.1/xm grains to velocities ~ 1 kms -1 (Yan et al. 2004). Grain-grain collisions 
at relative velocities ~ kms -1 will likely shatter high-porosity aggregates, and thereby limit the abundance of fragile 
high-porosity aggregates in the ISM. 

4 If Ay sa 10, the BAM2 model would be allowed by the Tx diagnostic if d/D < 0.24, but the TZx diagnostic requires d/D 0.50, so 
the BAM2 model would remain incompatible with the observational constraints. 
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